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Abstract 

The sensitivity of a pair of VIRGO interferometers to gravitational waves backgrounds of 
cosmological origin is analyzed for the cases of maximal and minimal overlap of the two 
detectors. The improvements in the detectability prospects of scale-invariant and non-scale- 
invariant logarithmic energy spectra of relic gravitons are discussed. 



1 Introduction and motivations 



Stochastic gravitational waves backgrounds |I] constitute a promising source for wideband 
interferometers whose operating window can be approximately located between few Hz and 
10 kHz. There are no compelling theoretical reasons why in such a frequency interval we 
should expect a negligible energy density stored in a stochastic background of primordial 
origin. Listening to phenomenology, we know that, unless Oqw (the logarithmic energy 
spectrum of relic gravitons) is either flat or decreasing (as predicted, for instance, by some 
classes of inflationary models), the present constraints on stochastic gravitational waves 
backgrounds are quite mild. Listening to the theory we know that if Oqw increases for 
frequencies larger than few mHz ||, then it is indeed possible to achieve a large signal in the 
operating window of interferometric detectors without conflicting neither with the fractional 
timing error of the millisecond pulsar's pulses || nor with the requirement that the total 
energy density of relic gravitons should be smaller than the total amount of relativistic matter 
at nucleosynthesis Q . These qualitative features can easily emerge in different models based 
on diverse physical frameworks including quintessential inflationary models 0, dimensional 
decoupling |§, early violations of the dominant energy conditions and superstring theories 

II- 

Recently the proposal of building in Europe an advanced interferometer of dimensions 
comparable with VIRGO jnj has been carefully scrutinized [11|. In view of this idea we 



would like to determine what would be the sensitivity of the correlation between two VIRGO- 
like detectors to a generic stochastic graviton background. The answer to this question de- 
pends upon two experimental informations (i.e. the specific form of the noise power spectra 
and the relative location and orientation of the two detectors) and upon the theoretical form 
of the logarithmic energy spectrum. Concerning the last point, the specific frequency depen- 
dence of ficw(/) is extremely relevant. In fact ficw(/) directly enters in the expression of 



the signal-to-noise ratio (SNR) J12|, 13] and, therefore, logarithmic energy spectra with a dif- 
ferent frequency dependence will lead, necessarily, to different SNR [f[4}| . Given a theoretical 
model, in order to compute reliably the sensitivity we have to specify the location of the two 
VIRGO detectors and the form of the noise power spectra. Thus, it is important to analyze 
how the relative distance between the two detectors of the pair can affect the sensitivity 
to the specific logarithmic energy spectrum of relic gravitons we ought to detect. In this 
respect, scale-invariant and non-scale- invariant logarithmic energy spectra lead to different 
sensitivity levels for the VIRGO pair already under the assumption that no improvement in 
the reduction of the thermal noises will take place prior to the construction of the second 
VIRGO detector. If a reduction in the contribution of the pendulum and pendulum's inter- 
nal modes to the noise power spectra is achieved the improvements in the sensitivity will be 
even sharper |15 . 
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2 Maximal and minimal overlap 



In the hypothesis that the graviton background is isotropic and unpolarized the reduction 
of the sensitivity due to the relative distance and orientation of the two detectors depends 
upon the frequency / and it is usually parameterized in terms of the (dimensionless) overlap 
reduction function |13| . 

By denoting (in spherical coordinates) with Q = (cos0 sin#, sin0 sin#, cos#) a generic 
direction along which the given gravitational wave (GW) propagates, the overlap reduc- 
tion function 7(/) can be expressed, in terms of the pattern functions determining the 
response of the i-th detector (i = 1, 2) to the A = +, x polarizations^ 



F A (f, n,i>) 



Tr{D(r) £ A (fi,^)} 



(1) 
(2) 



Notice that Ar = r\ — r 2 = ds is the separation vector between the two detector sites and the 
summation over the index A has to be taken over the physical polarizations (characterized 
by the polarization angle ip) 



e + (Q,t/j) = e + (fi)cos2?/> - e x (fi)sin2^ 
£ x (tt,ip) = e + (fi)sin2^ + e x (fi)cos2^ 



(3) 



of the incoming wave. The symmetric, trace- less tensor D(f) appearing in Eq. ([]]) depends 
on the geometry of the detector located in f. A general analytical expression for the function 
7(/) can be written as |13 



7 (/) = p (6) D{ 3 D 2ij + Pl {5) D{ 3 D 2 k t 8j s k + p 2 (S) D? D? s t Sj s k s t 
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with jk{S) the standard spherical Bessel functions: 
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1 We introduced the notation 
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to denote the average over the propagation direction 



and the polarization angle ip of the GW. 
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The normalization F is given by 



f = Y / (F 1 A (f 1 ,n^)F 2 A (f 2 ,n,^)^ | ls2 , (6) 

A 

where the notation 1 = 2 is a compact way to indicate that the detectors are coincident 
and coaligned and, if at least one of the two is an interferometer, the angle between its 
arms is equal to 7r/2 (L-shaped geometry). In this situation, by definition, j(f) = 1. When 
the detectors are shifted apart (so there is a phase shift between the signals in the two 
detectors), or rotated out of coalignment (so the detectors have different sensitivity to the 
same polarization) it turns out that |7(/)| < 1. The normalization factor F is 2/5 in the 
case of the VIRGO pair and in the case of any pair of wide band interferometers. 

By varying the frequency / the overlap reduction function is of order one until it reaches 
its first zero and, then, it starts oscillating around zero. The position of this first zero is 
roughly proportional to the inverse of the distance d between the two detectors. The location 
of the corner station of the first interferometer will be assumed in Cascina (43.6 N , 10.5 
E) where the VIRGO interferometer is presently under construction. The position of the 
second corner station is, at present, still under study. In our investigation we will suppose 
that the second corner station in different european sites. 

On a purely theoretical ground one ought to have a situation where the overlap reduction 
function is as close as possible to one for most of the frequencies in the operating window 
of the VIRGO pair. In other words we would like to push the first zero in the ultra-violet 
because this would imply that the region of maximal overlap (i.e. 7(/) ~ 1 ) gets larger. 
Since an increase in the region of maximal overlap produces a decrease in the relative distance 
of the two interferometers it will not be possible to decrease the distance ad libitum. In 
fact, when we decrease the distance between the detectors, we might introduce correlations 
between the local seismic and electromagnetic noises. We will assume that a distance of 
approximately 50 km is sufficient to decorrelate these noises. At the moment we do not have 
indications against such an assumption. 

The maximization of the overlap constitutes an important component as we can argue 
from the full expression of the signal-to-noise ratio. For the correlation of two interferometers, 
under the assumptions that the detector noises are Gaussian, much larger in amplitude than 
the gravitational strain and statistically independent on the strain itself, it can be shown 



12} |13| that the signal-to-noise ratio in a frequency range (f m , /m) is given, for an observation 
time T, by | 

3 ^o Pr\ f fM , f 7 2 (/)^w(/) V '' 2 



snr = ^p- vt / df 1 ; / uw ; 9 7 , (7) 

5v/2V \Jf m ps£\f)sP(f)) y> 

where H is the present value of the Hubble parameter. In Eq. (|7|), the performances 
achievable by the pair of detectors is controlled by the noise power spectra (NPS) S^' 2 \ 
whereas ficw(/) is the theoretical background signal defined through the logarithmic energy 

2 We follow the notations of Refs.^TJ, (TJ 
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spectrum, normalized to the critical density p c , and expressed at the present (conformal) 
time[] r/o 

Pc din/ 

It is intuitively clear, from the combined analysis of the two previous expressions, that if 
7(/) reaches its first zero a low frequencies, the value of the integral will be smaller than in 
the case where 7(/) reaches its first zero at larger frequencies. 



3 Overlap versus noise power spectra 



Given a specific logarithmic energy spectrum of primordial gravitons the signal-to-noise 
ratio is determined by the interplay between the overlap reduction function and the noise 
power spectra of the detectors. The noise power spectrum of the VIRGO detector can be 
approximated by an analytical fit [|l6|], namely 
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where 



10" 44 s 



fa 



500 Hz 



fb 



2 Hz 



E x = 3.46 x 10- 6 
E 2 = 6.60 x 10~ 2 
E a = 3.24 x 10~ 2 



The noise power spectrum of the VIRGO detectors is reported in Fig. |I[ In our parametriza- 
tion, Ei and E 2 control the contribution of the pendulum and pendulum's internal modes 
|T7|| to the thermal noise. E 3 controls instead the shot noise. For frequencies smaller than 



fb the noise power spectrum goes to infinity. The different values of Ej.,2,3 together with 
f a and fb define a specific noise configuration of the detectors. If, in the future, one of the 
contributions to the noises will be reduced, the noise power spectra will change accordingly 
and the noise configuration might be different. Without entering into the details of the 
actual experimental method which could allow such a reduction we can parametrize possible 
changes in the noise power spectra through a noise "reduction vector" 



P = (pi,P2,Ps) 

where pi < 1 defines the reduction in the corresponding coefficient Ej entering Eq. 
Within the present noise configuration p— (1, 1, 1). 



(10) 



3 In most of our equations we drop the dependence of spectral quantities upon the present time since all 
the quantities introduced in this paper are evaluated today. 
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Figure 1: We report the analytical fit of the rescaled noise power spectrum £„ defined in 
Eq. (|) in the case of the VIRGO detector. With the full (thick) line we denote the total 
NPS. We also illustrate the separated contribution of the (gaussian and stationary) sources 
of noise, namely the thermal (pendulum and pendulum's internal modes) and the shot noise. 



By using into Eq. (0) the expression of the theoretical spectrum given in Eq. (|8|) we have 
that the minimum normalization Q of a spectrum with functional variation oo(f), detectable 
by the VIRGO pair in an observation time T with a given SNR is determined by 

4.0 X10-' figy* , 



°" - J V T 
In this equation the information of the specific ui(f) is encoded in the quantity J 



where the integration variable is v — f /fo (/o is a generic frequency scale within the region 
fm 5: / < /m )■ We can assume / M = 10 kHz, whereas the lower extreme / m is put equal 
to the frequency /& entering Eq. (^). The choice of /o is purely conventional and in view of 
our discussion we took f = 100 Hz. 

Following the terminology of Sec. II the curve denoted by A in Fig. ^| corresponds to 
the maximal overlap. The minimal overlap is represented by the curve C where the location 
of the second VIRGO detector coincides with the present location of the GEO [18| detector. 
For completeness we also illustrate a possible intermediate overlap (profile B) corresponding 
the situation where the two detectors are roughly 500 km far apart. In principle, the effect of 
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Figure 2: We report the overlap reduction function(s) for the correlation of the VIRGO 
detector presently under construction in Cascina (43.6 N, 10.5 E) with a coaligned interfer- 
ometer whose (corner) station is located at: A) (43.2 N, 10.9 E), d = 58 km (Italy); B) 
(43.6 N, 4.5 E), d = 482.7 km (France); C) (52.3 N, 9.8 E), d = 958.2 km (Germany). The 
third site (C) corresponds to the present location of the GEO detector. 

a maximization (or reduction) of the overlap between the two detectors of the VIRGO pair 
is not independent on the analytical form of the logarithmic energy spectrum uj(f). Indeed, 
as we will show in the following Sections a reduction in the overlap has a mild effect if the 
logarithmic energy spectrum is scale- invariant. However, if the logarithmic energy spectrum 
is non-scale-invariant (and it increases with frequency) then the effect can be sizable. 

4 Scale-Invariant Energy Spectra 

Suppose, as a warm-up, that w(f) = 1 so that the logarithmic energy spectrum is strictly 
scale- invariant. Then according to Eq. ([11]) the sensitivity of the VIRGO pair can be 
computed. Let us firstly assume that the two detectors have minimal overlap (i.e curve C in 
Fig. 0) and let us suppose that the VIRGO detectors will have the NPS given specifically 
by Eq. (Q). In other words the reduction vector will have its present value p = (1, 1, 1). In 
this case, from Eq. fll2D we can determine the sensitivity of the VIRGO pair. By using now 
into Eq. (|TT|) the expression of the theoretical spectrum given in Eq. (§) we can determine 
the sensitivity of the VIRGO pair after one year of observation (i.e. T = 7r x 10 7 s) and with 
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SNR = 1 : 

^n~8.5xl(T 8 . (13) 

Suppose now to repeat the same estimate by assuming the maximal overlap. We want also 
to assume that the pendulum and pendulum's internal modes are suppressed by a factor of 
a hundred corresponding to a noise reduction vector p = (0.01,0.01, 1). Under this second 
set of assumptions we have that for T = 1 yr and SNR = 1 the sensitivity of the VIRGO 
pair becomes 

/ign~1.9x 10 -9 , (14) 

The two previous examples are quite extreme but a more complete analysis of the interplay 
between maximization of the overlap and noise reduction is illustrated in Tab. 1. In the first 



Table 1: The minimum detectable f2 for the three different location of the second VIRGO 
detector (see Fig. |2|) as a function of the reduction noise vector p. 
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column of Tab. 1 we report the different values taken by the reduction vector p whereas in 
the second, third and fourth columns (from the left) we indicate the sensitivities achieved 
by a VIRGO pair under the assumption that the second (coaligned) VIRGO interferometer 
is located, respectively, in the three positions specified by the three profiles A, B and C of 
Fig. 0. The first corner station is always assumed, in our estimates, to be in Cascina. 
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In spite of the fact that the example of an exactly flat spectrum is academic and oversim- 
plified f] we can draw, from our exercise, few interesting hints. We can see that if no noise 
reduction and minimal overlap are simultaneously assumed the sensitivity is, comparatively, 
smaller than in the case where a selective reduction of the thermal noises and maximal over- 
lap are postulated. If the reduction does not occur in the thermal components of the noise 
the improvement in the sensitivity is negligible, but still present. 

Our considerations were deduced only assuming a reduction in the thermal components 
of the noise. In principle, we should also consider the possible effect of a reduction in 
the shot noise. If the reduction is selectively applied to the shot noise, the improvement 
in the sensitivity is negligible This can be easily understood from the analysis of 

Fig. [I]: the shot noise starts being the dominant contribution to the NPS for / ~ 1 kHz, 
i.e. in a frequency region where the overlap begins to deteriorate (see Fig. H). As far as 
the achievable sensitivity levels are concerned, we can notice that the three location are 
practically indistinguishable. This statement is even more accurate if the reduction in the 
thermal noise components gets larger. 



5 Non- scale- invariant spectra 



It is interesting to repeat a similar analysis in the case where u)(f), instead of being scale- 
invariant, increases with /. In principle we would expect that in the latter case the impact of 
the maximization of the overlap will be more pronounced. The reason is that the contribution 
of uj(f) to the integrand of Eq. (|7|) (or of Eq. (|T2|)) increases at high frequencies if u(f) 
increases. Therefore, if the first zero of 7(/) falls just after 100 Hz (or possibly even before) 
the contribution of u>(f) will be erased more efficiently 

In order to show this behaviour let us analyze some specific examples among the ones 
mentioned in the introduction. For instance, in string cosmological models, the minimal 
pre-big-bang spectra have a two-branch form which can be expressed as 0, || 15 
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fx/ fs and g s are, respectively, the red-shift during, and the value of the 

r5J. The maximal amplified 



In this formula z s 

coupling constant at the beginning of, the string phase 



'Needless to say that the analysis of the Sachs- Wolfe contribution to the detected amount of anisotropy 



implies, on a theoretical ground that, for /„ 
originated during inflation. 
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if a scale-invariant spectrum was 
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frequency fi of the graviton spectrum is 

h (770) ~ 64.8 V^T ( — ] GHz, (17) 



10 3 X V12 



where n r is the effective number of spin degrees of freedom in thermal equilibrium at the 
end of the stringy phase, and g\ = M s /Mpi is the mismatch between the string (M s ) and 
Planck (Mpi) masses. The value of g\ corresponds to the dilaton coupling at the beginning 



of the radiation dominated epoch and it can be estimated to lie between 0.3 and 0.03 |L9 . 

In the case of increasing logarithmic energy spectra the evaluation of the sensitivity is a 
bit different from the case of purely flat spectra. If 00(f) grows in frequency the integrated 
spectrum can become, in principle large. We know, however, that the total amount of 
gravitons present at big-bang nucleosynthesis (BBN) cannot exceed the total amount of 
relativistic matter M. Otherwise the expansion rate of the Universe would increase too 
much and the observed light elements abundances could not be correctly reproduced. This 
implies that 

/"/max 

h 2 / n GW (/ J% ) din/ < ~ 5 x 10" 6 , (18) 



where f2 7 (^ ) = 2.6 x 10 is the fraction of critical energy density stored in radiation at 

the present observation time rj ; f ns ~ 10~ 10 Hz and / max are, respectively, the nucleosynthesis 
frequency and the maximal frequency of the graviton spectrum (for our present example 
/max — fi)- It is intuitively clear that if the spectrum is flat the BBN bound will be easily 
satisfied provided the theoretical amplitude of the spectrum, fl is roughly less than 10~ 6 . 
In the case of growing spectra the situation is more tricky. Let us denote with 



' max 



maY 5 x 10 6 ff m 

, 1= w(/)dln/, (19) 



1 1 h 

the maximal normalization of the spectrum compatible with the BBN bound. The sensitivity 
to a given 00(f) (which now increases with /) will be always given by Eq. (|TTj) . We should 
however exclude from the parameter space of our theoretical model those regions where 
VI > fi max . Thus, the regions of the parameter space of a given model for which f2 max > f2 
are, simultaneously, visible by the VIRGO pair and compatible with the BBN. In the absence 
of a specific theoretical normalization this would be the end of the story. However, in 
some models one can also estimate (with a given accuracy) the theoretical normalization 
of the spectrum. For instance, in the case of string cosmological models, the theoretical 
normalization of the spectrum can be expressed as 



Q th ~2.6gl\^-j 0,(770). (20) 

The accuracy of this determination coincides with the accuracy in the determination of g\ 
(the dilaton coupling at the beginning of the string phase) and of n r as defined in Eq. (|Ti 
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If we want to compare the achievable sensitivity not only with the BBN bound but also 
with the theoretical normalization we will have to require that f2 max > f2 and f2 > 
are simultaneously satisfied. With these necessary specifications let us now analyze the 
impact of the maximization of the overlap in the case of non-scale-invariant spectra of string 
cosmo logical type. Our results for the ratios Q max /Q and Q /Q, are reported in Figs. [| 
[|, and In Fig. [] we illustrate the sensitivity to string cosmological spectra under the 
assumption of maximal overlap (profile A of Fig. |2|). The region in the left plot of Fig. [3] 
traces the area of the parameter space where f2 max > £7, whereas the right plot corresponds 
to the region of the parameter space for which Q > Q. The plots of Fig. § can be compared 




Figure 3: We illustrate the regions of the parameter space of non-scale-invariant spectra 
of string cosmological type. As a function of g\/g s and z s we report the regions for which 
fi max /0 > 1 (left) and O /fl > 1. The overlap is assumed to be maximal (profile A of 
Fig. |2|). The fiducial set of parameters chosen for these plots (and determining the maximal 
frequency of the spectrum) is g\ = 1/20 and n r = 10 3 . 

with the case where the overlap is not maximal. This is done in Figs. |] and [|. In Fig. [| 
the meaning of the shaded regions is exactly the same as in Fig. |3| but with the assumption 
of minimal overlap. In fact, the relevant profile of 7(/) is the one labeled with C in Fig. 
0. By comparing Fig. [3] with Fig. ^ we can notice that the effects of the maximization 
of the overlap is to increase the visibility region both in terms of f2 max /f2 and in terms 
of Q /Q. As far as the maximal values of these two ratios are concerned the effect of the 
overlap reduction is irrelevant. Thence, we can conclude that a maximization (minimization) 
of the overlap between the VIRGO detectors is more evident if the detectable spectrum is 
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og z s log z s 

Figure 4: In order to make easier the comparison we report the same quantities discussed 
in Fig. fj, with the same fiducial choice of g\ and n T but under the assumption of minimal 
overlap, i.e. profile C in Fig. [2|. 
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non-scale- invariant. In the latter case, an increase in the overlap can sizably enlarge our 
sensitivity to the parameter space of a given model. 

For sake of completeness we also report in Fig. [| the patterns of the visibility region for 
non-scale-invariant spectra in the case of intermediate overlap. In this case the two VIRGO 
detectors are assumed to be roughly 500 km apart and this corresponds to the profile B of 
Fig. As we can see from Fig. |5] the area of the visibility region is, approximately, in 
between the ones obtained in the case of Figs. |3| and f|. 

In order to complete our analysis we would like to discuss the simultaneous effect of 
overlap maximization and noise reduction since this is one of the open possibilities of the 
upgraded VIRGO program. Suppose, for instance, that the noise configuration of the up- 
graded VIRGO would correspond to a reduction vector p — (0.1,0.01, 1). If the overlap is 
either maximal or minimal the visibility region are the ones reported in Fig. |]. In particular, 
in Fig. |^ we report the regions for which Q > Q (needless to say that we report only those 



n lh /n o. th /n 




6 8 10 12 14 16 IS 6 8 10 12 14 16 18 



og z s log z s 

Figure 6: Under the assumption of selective thermal noise reduction we illustrate the be- 
haviour of visibility region of the VIRGO pair in the case on growing logarithmic energy 
spectra of string cosmological type. We assumed a noise reduction vector p = (0.1,0.01, 1). 
At the left the overlap is maximal whereas at the right the overlap is minimal. As in the 
previous plots we consider a fiducial set of parameters (gi = 1/20 and n r = 10 3 ) to fix the 
maximal frequency. The sharp drop in the sensitivity after z s = 10 s can be understood by 
noticing that, in this case, / s ~ fo ~ 100 Hz. As we can see by comparing the left plot to 
the right one, if the thermal noises are consistently reduced, the gain in overlap gets less 
relevant. This is not the case if the noise are not reduced. 



12 



regions of parameter space for which the BBN bound is satisfied, i.e. Q th < fi max ). In the 
left plot of Fig. ^| we assumed maximal overlap (profile A of Fig. |2|) whereas in the right 
plot we assumed minimal overlap (profile C of Fig. |2|). As we can see, a joined reduction of 
thermal noises overwhelms almost completely the effect of the maximization of the overlap 
in the sense that, if the thermal noises are consistently reduced, the visibility region seems to 
be rather insensitive to the relative location of the two detectors. Similar conclusions (and 
similar plots) are obtained in the case of different reduction vectors affecting the thermal 
noises. Finally, if the shot noise contribution is reduced the sensitivity is only mildly affected 
[ 15 1 . The reason for this statement stems from the fact that the shot noise reduction starts 
being important for frequencies larger than the kHz (see Fig. |TJ) |15| . But for / > 1 kHz the 
overlap reduction is very efficient in spite of the location of the two detectors as it can be 
argued from Fig. 0. 



6 Conclusion 

The effect of maximization (or reduction) of the overlap of two coaligned VIRGO detectors 
depends upon the specific form of the logarithmic energy spectrum of relic gravitational 
waves backgrounds. If the spectra are purely scale-invariant the effect of minimization of the 
overlap are negligible. If the spectra are non-scale-invariant the maximization of the overlap 
has a sizable impact on the sensitivity. 

If we simultaneously maximize the overlap and reduce (selectively) the contributions of 
the pendulum and pendulum's internal modes to the noise power spectra, the sensitivity 
of the VIRGO pair to a scale-invariant spectrum can be as low as 10~ 9 after one year of 
observations and with SNR = 1. 

In the assumption that no noise reduction will take place in the context of the upgraded 
VIRGO program, the maximization of the overlap alone enlarges sizably the visibility window 
(always for T = 1 yr and SNR = 1) for the case of non-scale- invariant (and growing with 
frequency) logarithmic energy spectra. In the assumption of a consistent reduction of the 
thermal noises the sensitivity to growing logarithmic energy spectra increase significantly in 
spite of the location of the second VIRGO interferometer since, in the latter case, the good 
effect of the overlap maximization is overcome by the thermal noise reduction. 

From a theoretical perspective, our results support the conclusion that the maximization 
of the overlap has a different impact depending upon the analytical form of the logarithmic 
energy spectrum, urging, in particular dedicated studies of the forthcoming data for the 
specific case of growing logarithmic energy spectra. From a purely experimental perspective 
our findings suggest that a simultaneous maximization of the overlap and a noise reduction 
in the pendulum and pendulum's internal modes is desirable and extremely useful for a 
decisive improvement in the sensitivity of the VIRGO pair. 
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